Hamilton's turns as visual tool-kit for designing of single-qubit unitary gates 
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Unitary evolutions of a qubit are traditionally represented geometrically as rotations of the Bloch 
sphere, but the composition of such evolutions is handled algebraically through matrix multiplication 
[of SU(2) or SO(3) matrices]. Hamilton's construct, called turns, provides for handling the latter 
pictorially through the as addition of directed great circle arcs on the unit sphere S'^ C R^, resulting 
in a non-Abelian version of the parallelogram law of vector addition of the Euclidean translation 
group. This construct is developed into a visual tool- kit for handling the design of single-qubit 
unitary gates. As an application, it is shown, in the concrete case wherein the qubit is realized as 
polarization states of light, that all unitary gates can be realized conveniently through a universal 
gadget consisting of just two quarter- wave plates (QWP) and one half- wave plate (HWP). The 
analysis and results easily transcribe to other realizations of the qubit: The case of NMR is obtained 
by simply substituting 7r/2 and tt pulses respectively for QWPs and HWPs, the phases of the pulses 
playing the role of the orientation of fast axes of these plates. 

PACS numbers: 03.67.Lx, 03.65. Fd, 03.65.Vf, 42.50.Dv 
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I. INTRODUCTION 



States of a qubit are in one-to-one correspondence with 
points of the (solid) unit ball C {D for disk); de- 
pending on the context it is called the Poincare or Bloch 
ball. Pure states are on the boundary of and mixed 
states correspond to the interior points. The von Neu- 
mann entropy which is a measure of the mixedness of the 
state p has the simple form 

1 -I- r ^ -\- r 

S{P) = -Y. ^log2^, 0<r<l, 



where r is the radial distance of the point representing 
the state measured from the center of . When the 
qubit is realized as polarization states of hght, r equals 
the degree of polarization [l|, Q • The center corresponds 
to the maximally mixed (completely unpolarized) state. 

To go with this attractive geometric portraying of 
states, unitary evolutions p — >■ U pU\ U G SU(2) act 
as (three-dimensional) rotations leaving the center of the 
state space invariant. This is a realization of the ad- 
joint representation of SU(2) as the two-to-one SU(2) — > 
S0(3) homomorphism, both U and — C/ of SU(2) imag- 
ing to the same element of SO (3). More general physical 
evolutions or channels act as inhomogeneous linear maps 
on D^; in addition to mapping into itself, some ad- 
ditional requirement amounting to complete positivity Q 
will have to be satisfied by these maps. Such nonunitary 
evolutions, however, play no role in this work. 
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Though states and their (unitary) evolutions are thus 
represented geometrically, composition or concatena- 
tion of evolutions is traditionally handled algebraically 
through matrix multiplication. Hamilton's turns Q offer 
a visual tool for handling this last aspect too in a geomet- 
rical or vivid pictorial manner. In this picture, unitary 
evolutions are represented by (equivalence classes of) di- 
rected great circle arcs on S^ , with composition of unitary 
evolutions correctly represented by a geometric addition 
rule for these directed arcs, quite analogous to the man- 
ner in which translation group elements in an Euclidean 
space are composed using the parallelogram law of vector 
addition. 

An extensive description of Hamilton's construct may 
be found in the book of Biedenharn and Louck Q , while 
a simplified presentation of the addition rule for turns 
is given more recently in Ref. 0] ■ An early application 
of this construct to polarization optics can be found 
in , and generalization of the construct to other low- 
dimensional groups can be be found in . 

The principal aim of the present work is to develop 
Hamilton's geometric construct into a tool-kit for han- 
dling the composition and synthesis of unitary single- 
qubit gates in an efficient pictorial manner with no re- 
course to matrix multiplication. To be concrete, it is 
assumed in much of the presentation that our qubit is 
realized as polarization states of a photon |12h14| , but 
transliteration to other realizations of the qubit will be 
evident. For instance, the case of NMR quantum compu- 
tation is easily seen to correspond to quarter- wave plates 
(QWPs) and half-wave plates (HWPs) being replaced, re- 
spectively, by 7r/2 and tt pulses, the phases of the pulses 
playing the role of the orientation of the fast axis of the 
plates. 

The tool-kit presented has the obvious limitation that 
it can handle only single-qubit (unitary) gates. Even so 
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it could prove to be of value to quantum computation in 
view of the fact that single-qubit gates, along with just 
one two-qubit gate like the C-NOT gate, can realize all 
multi-qubit gates [Tsl - fisj . 

The usefulness of this tool-kit is not limited to the 
domain of quantum information and computation. The 
group SU(2) pervades many areas of science, either di- 
rectly or through the rotation group S0(3) — SU(2)/Z2, 
and this tool-kit could therefore prove useful in these 
other areas as well. In particular, it is of direct interest 
to classical polarization optics. 

The entire presentation shoots towards the main result 
formulated as a theorem at the end of the paper, which 
asserts that all single-qubit gates can be conveniently 
realized using a universal gadget consisting of just two 
QWPs and one HWP. We should hasten to add, however, 
that this theorem is not new in itself, but has been formu- 
lated earlier using algebraic methods [Tol . and Bagini et 
al. (20j have presented a particularly helpful exposition of 
this result of Ref. [l^ which has been variously used [2l| - 
[27| . Whereas the algebraic approach took a sequence of 
several papers 0, HM l29[ to eventually arrive at the fi- 
nal result in the fourth [1^ , through a sequesnce of false 
starts and refinements, the pictorial approach presented 
here will be seen to render the result visual, and almost 
obvious. 

Since this geometric construct of Hamilton, called 
turns does not appear to be as well known as it deserves 
to be, we begin with a description of this construct it- 
self, relating it to the three prominent parametrizations 
of SU(2) — the homogeneous Euler, the axis-angle, and 
the Euler parametrizations — and bringing out its inter- 
esting connection with the Berry-Pancharatnam geomet- 
ric phase. 

II. HAMILTON'S TURNS 

Reversible gates acting on a qubit are in one-to-one 
correspondence with 2x2 unitary matrices u £ SU(2). 
The SU(2) matrices can be conveniently described by any 
triplet Ti,T2, T3 of Pauli-like Hermitian matrices satisfy- 
ing the defining algebraic relations 

TkTl — 6klTo + ieklmTm, (1) 

where tq = l2x2 is the unit matrix £ SU(2). To be 
specific, we take these matrices to be ri = cr3,T2 — a\ 
and T3 — where (Tj's are the standard Pauli matrices. 
The family of all unitary (reversible single-qubit) gates 
u € SU(2) then get parametrized as 

(ao — ia\ — m2 — as \ 

Oq + + 03 + flg = Aq + a • a = 1. (2) 

That is, the four real parameters (cq, a) = (ao, ai, 02, 03), 
called the homogeneous Euler parameters, correspond to 



a point on the three-sphere S'^ C M^. Elements of SU(2) 
are thus in one-to-one correspondence with points on S'^, 
consistent with and exhibiting the fact that S'^ is the 
group manifold of SU(2). Hamilton's turns constitute a 
powerful visual representation of this C on C 
through (equivalence classes of) directed great circle arcs, 
with group multiplication of SU(2) matrices (concatena- 
tion of single-qubit unitary gates) faithfully transcribed 
into a "parallelogram law of addition" for these directed 
geodesic arcs on S^. 

Given a unitary gate u or (ao,a) € S'^, the constraint 
a^ + a- a = \ guarantees that wc can find an ordered pair 
of unit vectors ni, 712 G M'^ such that 

ao = ni • n2, a = niAn2- (3) 

It follows that a directed great circle (or geodesic) arc 
on S^, with tail at fii and head at rj.2, can be associated 
with the unitary gate u. Clearly, such an association is 
not unique, and the non-uniqueness is precisely to the 
following extent: Any pair h'l, n'2 obtained by rotating 
both fii and 7x2 by equal amount about ni An2 will meet 
the requirements in Eq. ([3]) , and hence will correspond to 
the gate represented by the original pair hi,n2- Such 
a rotation obviously corresponds to rigidly sliding the 
directed arc representing u along its great circle. 

One is thus led to consider equivalence classes of di- 
rected great circle arcs on S^ , the equivalence being with 
respect to the sliding just noted: Two such directed arcs 
are equivalent if they are on the same great circle and 
if one can be made to coincide with the other by rigidly 
sliding it on the great circle. These equivalence classes 
are called Hamilton's turns. It is clear that elements of 
SU(2) are in one-to-one correspondence with turns (as- 
suming the arclength of turns is restricted not to exceed 
^). 

If u € SU(2) is represented by the turn whose rep- 
resentative element is the directed great circle arc from 
ni(tail) to n2(hcad), we denote this fact through 

u ~ T{hi, 77,2) = ni • 712 — ifii A 77,2 • r. (4) 

Henceforth wc talk of turns. SU(2) matrices, and unitary 
single-qubit gates interchangeably. For brevity, we often 
call a representative arc itself as the turn and its length 
as the length of the turn or simply as the turn length, but 
this abuse of terminology should cause no confusion. 

Two special elements of SU(2) namely tq and —tq are 
distinguished in that they constitute the center of the 
group. This distinction should be expected to manifest 
itself in any representation, and the one due to Hamilton 
happens to be no exception. The unit clement tq, the 
trivial gate, corresponds to the null turn ni = 77,2 G 
S^, and the gate —tq to ni —77,2 G S^, the respective 
turn-lengths being 0, tt. The equivalence class associated 
with either is clearly a two-parameter family, since hi 
can be any point on S^. However, the equivalence class 
of directed great circle arcs associated with any other 
turn is a one-parameter family. In particular, every turn 
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or It € SU(2), u ±To, has associated with it a unique 
directed great circle of . 

Analogy with the Euchdean translation group (in two 
dimensions, for instance) wherein group elements are rep- 
resented by equivalence classes of ^free vectors^ is obvious. 
There the abelian group composition takes the geometric 
form of parallelogram law of vector addition. It turns out 
that such a geometric or pictorial composition of group 
elements applies to the present non- Abelian case of SU(2) 
as well, with turns playing the role of (equivalence classes 
of) free vectors; indeed, the power of Hamilton's turns 
can be traced to this fascinating fact. 




FIG. 1: The "parallelogram law" or "addition rule" for turns, 
with turn AB = turn BE and turn BC = turn DB representing, 
respectively, unitary gates u and v and turn AC and turn DE repre- 
senting, respectively, the products vu and uv. Since turn LK repre- 
sents vu and turn KM, the inverse of turn DE, represents v~^u~^, 
and turn LM represents the commutator [v, u] = v~^u~^vu. 

To see this pictorial composition law, assume that we 
are given two SU(2) gates u, v and we wish to compute 
geometrically (visually) the matrix product vu. Referring 
to Fig. [3 let the directed great circle arc AB represent 
u and let BC represent v. It is important to note that 
the representative arcs are so chosen that the head of the 
right factor u and the tail of the left factor v coincide at 
B^ since great circles on certainly intersect, this 
can always be arranged for any given pair of turns. Now 
draw the directed geodesic arc from the free tail A^ ni 
to the free head C ~ to obtain a new turn repre- 
sented by the directed arc AC. The important claim is 
that turn AC correctly represents the matrix product vu. 
This is readily verified through matrix multiplication: 

vu = T{n2,n3)T{ni,h2) 

= (n2 • ns — in2 A ■ T)(ni • n2 — irJ-i A n2 • r) 
= ni • fia — mi A n3 • X = T(ni, ng). (5) 

Remark: The only property of the r-matrices used in 
this verification is TkTi = SuiTq + iekimTm, and hence this 
result and its consequences are independent of which set 
of Pauli-like matrices was actually used. 



We forced ourselves to perform the kind of matrix mul- 
tiplication in Eq. ([5]), just once, simply to demonstrate 
that the one-to-one correspondence between SU(2) gates 
and turns is indeed a group isomorphism. The rest of 
this work, however, will rest solely on the geometric or 
visual construct of turns, with almost no recourse to ma- 
trix multiplication. We note in passing that the above 
composition law immediately implies that the matrix in- 
verse of the SU(2) gate or turn T{hi,n2) corresponds to 
r(n2,ni), the reversed turn. 

To compute or construct the 'other' product uv, choose 
E and D on the great circles, respectively, through A, 
B and C, B such that AB = BE and BC = DB. Then 
turn BE = turn AB represents u and turn DB = turn BC 
represents v. It is thus obvious in view of Eq. ([5]), that 
turn DE represents the product uv, which is manifestly 
different from the product vu (represented by turn AC), 
giving a vivid pictorial depiction of the noncommutative 
nature of the "addition rule" for turns, consistent with 
the non-Abelian nature of SU(2) composition. 

Remark: In spite of being quite different from one an- 
other, turn AB and turn DE share one important com- 
mon aspect. To see this, consider the spherical triangles 
ABC and EBD. The angle at B is the same for both tri- 
angles, AB = BE, and CB = BD. The two triangles are 
thus congruent, showing that turn AC and turn DE have 
the same turn length. As we shall see, this is a pictorial 
manifestation of the fact tr vu = tr uv. 

Presented also in Fig. [T] is a visual display of the com- 
mutator of two SU(2) gates. Recall that the commu- 
tator of a pair of elements gi, g2 of a multiplicative 
group G is defined as [g2, gi] = 92^9i^929i, the mul- 
tiplicative difference of (71(72 and (7251- Let the geodesic 
arcs DE and AC when extended meet at K. Choose 
points L, M on these extended arcs such that AC = LK 
and ED = KM. Since turn EB corresponds to u^^ and 
turn BD to turn ED corresponds to the product 

v~^u~'^ = {uv)~^. Referring now to the spherical triangle 
LKM, since turn LK corresponds to vu and turn KM to 
v~^u~^ we deduce, again in view of Eq. (O, that turn LM 
corresponds to the product v~^u~^vu, the commutator 
[v, u] of interest. 

It is often convenient to rewrite the SU(2) composition 
rule of matrix multiplication 

T(n3,n4)T(n2,n3)r(ni,n2) =T(ni,n4) (6) 

as the geometric (visual) rule of " addition" of turns 

turn hih2 + turn n2i^3 + turn h'ifii = turn fiifi^. (7) 

In transcribing from the " multiplication" mode of Eq. (|6]) 
to the "addition" mode of Eq. ([7]), however, it is impor- 
tant to remember that the individual terms in this " sum" 
in Eq. ([7]) read from left to right correspond to the fac- 
tors in the SU(2) matrix product Eq. ([6]) read from right 
to left] the order is important, the "sum" being non- 
commutative. This geometric "addition rule" for turns, 
which is clearly reminiscent of the parallelogram law for 
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the composition of elements of the (Abelian) Euchdean 
translation group, is associative and faithfully represents 
the non- Abelian or non-commutative group composition 
in SU(2). 

Our consideration of turns so far has been based on 
the homogeneous Euler parametrization of SU(2) Eq. (2). 
The group SU(2) can also be parametrized in the axis- 
angic form 



u = u{n, a) 



u{n, 2-k) 
u{h, 2tt + a) 



exp(— z— n • t) 

cos — Tn — i sm — n ■ T\ 
2 2 

-To, ii(n,47r) To, 
u{—n, 27r — a), 



(8) 



where n is a unit vector G M^. We may, in view of the 
last line of Eq. (|8]), restrict a to the range < a < 27r. A 
view of turns which corresponds to this parametrization 
proves more convenient for some purposes, and so we 
describe it briefly. 

We have seen that every turn, other than the special 
turns r(n, n) and T{n, ~n) corresponding to elements 
in the center oi SU(2), has associated with it a unique di- 
rected great circle; this great circle and an angle (length of 
the representative arc) fully specifies the turn. However, 
directed great circles and directed axes (or unit vectors 
n £ S^) are in one-to-one correspondence: if the directed 
great circle is specified by the ordered pair (ni ,712) of lin- 
early independent unit vectors, the directed axis is speci- 
fied by the unit vector n e 5^ in the direction of ni An2. 
We may thus denote a turn alternatively by the symbol 
r(n, a), with the understanding that the representative 
directed arc is on the great circle orthogonal to n and 
has arclength a. In other words, T(n, a) corresponds to 
u{n,2a) e SU(2). With the two special turns excluded, 
it is clear that this representation is unique [hi A n2 is 
nonvanishing for every u not in the center of SU(2), that 
is, for every turn whose turn- length £ (0, tt) ]. 

We note from the last line of Eq. (|8]) that T{n, TT + a) = 
T(— n, TT — a), and hence the restriction of a, the turn- 
length, to the range < a < tt. For the two special 
turns corresponding to the center of SU(2), we see that 
r(n, 0) and r(n, tt) represent, respectively. To and — to, 
independent of n G . 

Wc note in passing a direct relationship between the 
turn length and the trace of the associated SU(2) matrix. 
Indeed, the facts that tr w(n, a) = 2 cos ^ and u(fi, 2a) G 
SU(2) is represented by the turn T(n, a) of length = a 
show that tvu = 2cos£„. In particular, tru is positive 
or negative depending on whether the turn length ^„ is 
< or > 7r/2 : All traceless SU(2) matrices correspond 
to turn- length 7r/2, the popular Hadamard gate [l5| and 
HWPs considered below being examples. 



III. REALIZATION OF QUBIT AS 
POLARIZATION STATES OF LIGHT 

As indicated earlier, our illustrations demonstrating 
the power of turns in solving problems of synthesis of 
unitary gates will use the concrete context wherein qubit 
is realized as the polarization states of light. However, 
it is evident from the treatment to follow that the en- 
tire analysis applies equally well to other realizations of 
qubits and SU(2) gates, like NMR quantum computa- 
tion [15| or passive (lossless) linear optics of a pair of 
radiation modes at lossless beam splitters (s^l ■ 

Birefringent media play a particularly dominant role 
in polarization optics, both classical and quantum. Con- 
sider a (quasi-monochromatic) light beam propagating 
along the positive X3 direction of a Cartesian system 
{xi,X2,X3). The components Ei,E2 of the transverse 
electric field along the xi,X2 directions can be arranged 
into a column vector E G C^, called the Jones vector 
of the polarization state of the beam [l| . The intensity 
equals l-Bip -I- = E'^E. A Hncar optical system is 

correspondingly represented by a 2 x 2 numerical matrix 
J called the Jones matrix, and the input-output relation- 
ship is represented by 

^'in —J' £^out = JEin- (9) 

Lossless linear systems conserve intensity: E^^^Eo^t = 
El^Ein- It follows that the Jones matrices of such sys- 
tems are unitary. Birefringent media, which introduce 
a relative phase between a characteristic pair of orthog- 
onal linear polarization states, and optically active me- 
dia which introduce a relative phase between the two 
(orthogonal) circular polarization states are examples of 
such lossless linear systems of interest to polarization op- 
tics. Suppressing an overall phase, the Jones matrices of 
lossless linear systems can be identified with elements of 
the unimodular unitary group SU(2). In the case wherein 
the qubit corresponds to the polarization states of a pho- 
ton, these arc indeed the relevant unitary or reversible 
single qubit gates. 

A birefringent plate (compensator) whose "fast axis" 
is along the transverse xi direction has the Jones matrix 



J = Coiv) 



e*''/^ 



exp(-z|Ti) eSU(2), 



(10) 



77 being the relative phase introduced by the plate; we 
have rj ~ el/\, where £ is the thickness of the plate, e > 
is the difference between the refractive indices for the 
two characteristic orthogonal linear polarizations, and A 
is the wavelength. It is clear that if the fast axis is at an 
angle ip with the xi axis, then the Jones matrix would 
be 



~ cos{ri/2) To — i s'm{rj/2) (cos(2(^) ti 



■ sm{2ip) T2), 
(11) 
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where the two-dimensional matrix 

4>(^)=expH^r3)=f^°^^ '^"^^"j (12) 

y o, y smip COS If J ' 

is an element of the subgroup S0(2) C SU(2). 




FIG. 2: The sphere of turns T as it applies to the polarization 
qubit. All turns on the equator (horizontal turns) correspond to op- 
tical rotators, and vertical turns correspond to bircfringent plates. 
While relating to the subscripts of HWPs and QWPs it should be 
remembered that spatial rotation of magnitude ip/2 is reflected on 
T as polar rotation of magnitude ip and not ip/2. 

Quarter-wave plates (QWPs) and half-wave plates 
(HWPs) are particular cases of birefringent plates and 
correspond, respectively, to 77 = tt/2 and tt; they could 
equally well be called A/4 and A/2 plates. The Jones 
matrix of a QWP with fast axis along the xi -direction 
is thus cxp(— i^Ti) = (to — iri)/-\/2; we denote this 
QWP by Qo, so that = $(<p)Qo*(<^)"^ = [tq - 
i{TiC0s2ip + T2 sin 2(^) ]/-\/2 represents the QWP whose 
fast axis makes an angle ip with the positive xi axis. Sim- 
ilarly, we use the notation Hq for the HWP —ir^ so that 
stands for <i>((y9)iJo<I'(v)~^ = — i(Ti cos2iy9-f T2 sin2(y9), 
a HWP whose fast axis makes angle (p with the xi axis. In 
particular, H^^^ corresponds to the Hadamard gate . 
Thus H^p = [QipY, for all ip. These arc particular cases 
of a more general and evident fact: If an optical system 
represented by Jones matrix J is physically rotated by an 
angle about the positive x^ axis, the resulting system 
will have Jones matrix exp(— it/sra) Jexp(i93T3). 

The S0(2) matrix exp(— zt/JTs) plays yet another role 
in polarization optics: It is also the Jones matrix of an 
optically active medium. Specifically, an optically active 
medium (or simply optical rotator) which introduces a 
relative phase a between the left and the right circularly 



polarized states has the Jones matrix 

i?(Q;) = exp(-i-T3). (13) 

Numerically, R{a) = $(a/2); however, we have chosen to 
use a different symbol R{a) for the optical rotator to dis- 
tinguish it from $(•), which stands for physical rotation 
of a gadget in the transverse plane. 

We depict in Fig. [2] the sphere of turns T. It is clear 
that all "vertical turns" correspond to birefringent plates: 
QWPs and HWPs have turn lengths of 7r/4 and 7r/2, 
respectively. Turns on the equator correspond to optical 
rotators. 

In addition to the homogeneous Euler and axis-angle 
parametrizations of SU(2) gates already considered there 
exists a third one, the Euler angle parametrization, 

= exp(-i-^r3) exp(-z-77ri) exp(-i-CT3). (14) 

This paramctrization of the SU(2) group of unitary gates 
can be viewed as saying that every such gate is equiv- 
alent to an appropriate birefringent plate Co{r/) sand- 
wiched between two appropriate optically active media 
R{^), R{C), with the effective thicknesses of the three me- 
dia engineered to match the Euler parameters ^, 7?, C of 
the gate under consideration. However, from an experi- 
menter's point of view this cannot be the most convenient 
realization of the various unitary gates J € SU(2); unlike 
the HWP and QWP, Co (77) is not a component readily 
available off the shelf, and a variable R{£,) tends to in- 
troduce ^-dependent loss and dynamical phase. It turns 
out that QWPs and HWPs alone are sufficient. As we 
will see, the geometric representation of Hamilton ren- 
ders this fact particularly transparent and visual. 

We may note in passing that the Euler angle 
paramctrization [Eq. (jl4p ] can be rewritten in the modi- 
fied form 

«(e, '7, C) - c^niv) m + = m + c^cMv)- (15) 

This means an arbitrary unitary gate is a variable bire- 
fringent plate preceded or followed by a variable rotator. 
The fact still remains that, unlike QWPs and HWPs, 
variable birefringent plates are nonstandard polarization 
optical components, and optically active media introduce 
undesirable losses and dynamical phases which vary with 
the optical rotation or effective thickness of the medium. 

IV. HAMILTON'S TURNS AND 
TRANSFORMATION OF POLARIZATION 
STATES 

Polarization states are conveniently described (even in 
the classical case) by restricting attention to Jones vec- 
tors of unit norm (unit intensity) and ignoring an overall 
phase. Such normalized Jones vectors correspond, in the 
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quantum case, to state vectors of a two-level system or 
qubit. In either case, the state represented by a normal- 
ized Jones vector E is fuhy determined by the (complex) 
ratio z = E2/E1. This is rendered particularly transpar- 
ent by going over to the coherency or density matrix, and 
one then finds that polarization states are in one-to-one 
correspondence with points on the unit sphere , called 
the Poincare or Ricmann sphere, obtained by identifying 
the points z — 00 (the one-point compactification) of the 
complex z plane: 

ti^EE'') = 1 and e"^E{m) - E{rh), rh e S'^ 
p{rh) EE E{m)E{rh)\ 

Written in more detail. 



E{m) = 



1 



cos I + e 



V2 



§-e-/^sin§] 



(16) 



so that for any state m e we have for the ratio z{ifi) 
of the components of E{rh) the expression 



E. 



sm a sm ip + i cos f 



El 1 + sin 9 cos ip 
The corresponding coherency or density matrix reads 
p{m) = E{rh)E{rhy 



(17) 



1 + sin 9 cos (fi sin 9 sin (f — i cos ( 
sin 9 sin ip + i cos 9 1 — sin 9 cos tp 



z{m) = pi2(m)/pii(m). 



(18) 



The parameters 0, ip are respectively the 
polar and azimuthal coordinates of m = 
(sin 6* cos sin sin cos6') g S*^. 

Had we used in place of the " r matrices" the standard 
Pauli a matrices, the coherency matrix would have read 



PVm) = -{ao+m-cr) = - 



1 



smf 



1 — cos ( 



[ctq again is the unit matrix l2x2 = tq.] 
Eirh) would have been parametrized as 



E{rh) 



cos 9/2 



sin 6*72 



(19) 
Consequently, 

(20) 



so that the ratio E2/E1 becomes z{rh) = e^'^ tan(6'/2), 
a form more familiar in the context of NMR and the 
associated Bloch sphere. 




LCP 



FIG. 3: The polarization states of light represented as points on 
the Poincare sphere V. Right and left circularly polarized states 
occupy the polar positions, and points on the equator correspond 
to linear polarization states. Spatial rotation of the plane of polar- 
ization by angle ip corresponds to 2(/3 rotation about the polar axis. 
Thus linear polarization LPq along the xi axis and linear polariza- 
tion at 45° to the x\ axis are separated on V by 90°. Points 
of the upper (lower) hemisphere correspond, respectively, to right 
(left) handed elliptic polarization. 



Staying with the choice r, rather than a, we sketch 
in Fig. |3] the Poincare sphere V . We adopt the conven- 
tion that right and left circular polarization (RCP and 

LCP) states correspond, respectively, to ^ 



1 

V2 



and 



1 

72 



1 



or, equivalently, to the Poincare sphere co- 



ordinates (0,0,1) for RCP and (0,0,-1) for LCP. Lin- 
ear polarization at angle a to the (positive) Xi axis has 



Jones vector 



cos a 
sin a 



and corresponds to the Poincare 



sphere coordinates (cos 2a, sin 2a, 0). Thus, all linear 
polarization states are on the equator; all states above 
the equator have right-handed elliptic polarization, while 
those in the lower hemisphere are left-handed. Antipodal 
points of V correspond to orthogonal polarization states: 
E{—rh)'^E{rh) = 0, equivalently tr ( p(— m)/9(m) ) = 0. 

Remark on convention: Since the expressions for 
E{rh) and p{rh) may "appear to be" considerably sim- 
pler with the choice cr, one may wonder why one chose 
T in the first place. The reason is one of convention, a 
price one occasionally pays for tradition. In the case of a 
spin-i particle of the NMR context or a two- level atom. 



one chooses the vectors 



and 



to correspond 
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to the poles of and prefers to associate Bloch's name rotations: 
with this sphere. Obviously, these so-called computa- 
tional basis states are eigenstates of 0-3. In polarization t ■ 

optics, on the other hand, the circularly polarized states ^^^\ ~* '''1 ) ^ 72-1 (ai) - 

( ±1 ) given the special honor of polar positions, 

1 ^1 ^ • r n 11 • 1- • / cosq;2 sina2 

but these are eigenstates of cr? = 73. _bollowmg tradition 1 ■ '^'^ \ t-, / \ \ n 1 n 

,1 , • 1 • 1 , , • 1 expi — I — T2 ) — > 7C2 02 = 10 

we wish to keep the polar axis as the vertical ana third 2 \ — sin a cos a 

axis for polarization qubit. Cyclic permutation of the a ^ 

matrices seems to be the minimal way of meeting these / q,^ _ gjj^ q,^ q 

concerns or requirements without offending in any way gxp( —i — Ts)—^- 7^,3(03) = ( sina3 cosq!3 

the basic algebra [commutation and anti-commutation 2 \ 1 
relations, Eq. ([T]) ]. 



(23) 



ACTION OF TURNS ON THE POINCARE 
SPHERE AND CONNECTION WITH 
GEOMETRIC PHASE 



The notation T(n, a) for turns, corresponding to the 
axis-angle parametrization u{h, 2a) = cxp(— ian • r) of 
SU(2), proves convenient in exhibiting the action of turns 
on the Poincare sphere, our state space. Unitary gates 
T{n, a) act on the Poincare sphere, that is, on state 
p{fh) = ^[tq + rh ■ t] represented hy rh £ V, in this 
manner: 

T{n,a): p{rh) ^ p(m') = T(n, a) p(m) T(n, a) 

= exp(— i ah ■ t) p{rh) exp(i a n • x). 



(21) 




As expected, the two special gates tq, — tq corresponding 
to elements of the center of SU(2) are seen to have no 
effect on the Poincare sphere and, as noted earlier, the 
set of parameters (n, a) in T(n,a), with < a < tt. is 
unique for every other turn. 

Now, in view of the algebraic properties of the r matri- 
ces [Eq. dl])], the above transformation law simply reduces 
to 

T{n,a) : m ^ m' = {m ■ n)h + cos 2a [m — (rh • n)n] 

-I- sin 2a n Am. (22) 

Note that {rh • n)n is the component of rh along h and 
rh~ {rh-n)n is the component of rh orthogonal to n and 
hence in the plane spanned by rh, fi. Further h/\rh which 
is orthogonal to both rh and h has the same magnitude 
as m — [rh ■ n)n. Thus, the effect of T{n, a) on the 
state space or Poincare sphere V is an SO (3) rotation 
about fi, of extent 2a = twice the length of the turn. In 
particular, the three one-parameter subgroups of SU(2) 
gates exp( —i^Tk), fc = 1, 2, 3 , act on V, respectively, 
through the following one-parameter subgroups of S0(3) 



FIG. 4: Shown on the sphere of turns T is a geodesic tri- 
angle corresponding to a unitary resolution of identity. While 
turns AB, BC, CA compose to unity, their respective square roots 
compose to turn B'C whose turn length is half the area of the tri- 
angle. 

Now consider the spherical triangle ABC on the sphere 
of turns T shown in Fig.|4l the coordinates of A,B,C be- 
ing respectively (0, 0, 1), (1,0,0), (0,1,0). Since turn AB 
corresponds to the unitary gate —1x2, turn BC to —1x3, 
and turn CA to — iri, the (closed) triangular circuit, 

turn AB + turn BC + turn CA = null turn (24) 

represents a visual depiction of the matrix product iden- 
tity 

(-ZTl) {-iT3) i-lT^) = To (25) 

involving the Pauli gates — , a unitary resolution of the 
identity. 

Let A', B', and C be the midpoints of AB, BC, and CA 
respectively. Thus, turn AA' = turn A'B, turn BB' = 
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turn B'C and turn CC = turn C'A are the square roots 
of turn AB, turn EC, and turn CA, respectively. It is 
instructive to compute the composition of these square- 
root turns: turn AA' + turn BB' + turn CC. Note that 
AA', BB', and CC have equal arclength of 7r/4, and 
correspond respectively to -^{tq — iT2), -^(^0 ~ i^s), 
and -^{to ~ JTi). Since turn AA' = turn A'B, wc 
have turn AA' + turn BB' = turn A'B'. To compose 
turn A'B' with turn CC extend A'B' to meet the great 
circle through A and C at D. Comparing the spheri- 
cal triangles BB'A' and CB'D, we see that BB' = B'C, 
angle A'BB' = angle B'CD tt/2), and angle BB'A' = 
angle CB'D. The two triangles are thus congruent, and 
so turn A'B' = turn B'D and arclength of DC = ar- 
clength of A'B 7r/4). Thus, turn A'B' + turn CC = 
turn B'D + turn DC = turn B'C = turn BB'. Wc have 
thus shown 

turn AA' + turn BB' + turn CC = turn BB', 
i.e., y/-iTi \/-iT2 = \/-iT3 = -^(ro-ira). 

(26) 

The square roots thus compose to produce neither the 
null turn nor its square root (any turn of turn length 
=7r), but "^(''o ^ ^''"3)1 8- primitive eighth root of the null 
or identity turn tq. 

This is due to the following fact: Unlike triangles in 
the Euclidean plane, there is no notion of similar trian- 
gles (beyond congruence) in the spherical case; the area 
of a spherical triangle is fully determined by its three 
angles. This situation described by Eqs. (|25)) and (|26)) 
should be contrasted with the corresponding situation 
in respect of the Euclidean parallelogram law, wherein 
if three elements of the translation group compose to 
produce the null (identity), element then their respec- 
tive "square roots" too will certainly compose to yield 
the null clement, corresponding to a similar triangle with 
one-fourth area and the same interior angles as the orig- 
inal one. The failure of the SU(2) turns in this respect, 
as depicted by Eq. (PE)) . is rooted in the non-Abelian na- 
ture of the group on the one hand and in the nontrivial 
curvature of (as compared to the Euclidean plane) on 
the other; indeed, these two aspects go hand in hand, 
as may be seen also by consideration of the geometric or 
Pancharatnam phase. 

Before we turn to the geometric phase, we note, how- 
ever, that the sum of the three (interior) angles of the 
triangle ABC is in excess of tt by 7r/2, and this spherical 
excess equals the area of the triangle. We note also that 
the turn length of the composite turn AA' + turn BB' + 
turn CC, which is clearly a measure of the extent to 
which the square roots fail to compose to the null turn, 
is 7r/4, precisely half of the area of the original (closed) 
triangle ABC. That is, in the axis-angle notation T{n, a), 
the value of a corresponding to composition of the square 
roots (that is, turn B'C) is 7r/4; and n corresponds pre- 
cisely to fiA, the "starting point" A of the triangle. 



All these aspects apply not only to the particular tri- 
angle shown in Fig.|3]but also to an arbitrary geodesic tri- 
angle; construction of proof in the general case is slightly 
more elaborate, but very similar to the one for the spe- 
cial triangle in Fig. |4] (see Ref. [3l|). That is, we have for 
any spherical triangle ABC with area A 

turn AA' + turn BB' + turn CC = T{nA, A/2), (27) 

where A', B', C are, respectively, the midpoints of AB, 
BC, CA and fiA is the unit vector pointing in the direc- 
tion of A e 5^. 




FIG. 5: Generalization of the area formula of Fig.|4]to the case of 
a geodesic quadrilateral. 

A simple argument implies that this area formula ap- 
plies indeed to all geodesic polygons. Consider, for in- 
stance, the (geodesic) quadrilateral ABCD on T shown 
in Fig.m Let A', B', C, D' be the mid points of AB, BC, 
CD, and DA, respectively, and let E' be the midpoint of 
the geodesic CA. Let Ai be the area of the triangle ABC 
and A2 that of ACD so that A = A1+A2 is the area of the 
quadrilateral, and let be the unit vector correspond- 
ing to the point AeV. Wc have the closed-circuit rela- 
tion turn ABH-turn BC-fturn CD+turn DA = null turn. 
We wish to compose the square roots of these turns in 
that order. The basic idea is to break this quadrilateral 
into two triangles ABC, ACD (making use of the fact 
that turn CE' -I- turn AE' = null turn). We have 

turn AA' + turn BB' + turn CC + turn DD' 
= (turn AA' + turn BB' + turn CE') 

+ (turn AE' + turn CC + turn DD') 
= T{hA,Ai/2) + T{nA,A2/2) 
= T{nA, AI2). (28) 

In the last step we used the fact that all turns r(n, a) 
having the same ft form an (Abelian) one-parameter sub- 
group: 

r(n, a)T(n, a') = T(n, a')T{n, a) = T{n, a-^a'). 

(29) 
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Generalization to n-sided geodesic polygons is obvious. 

Remark: As a subtle but important aspect we may note 
that the addition rule in Eq. ()29|) and its covariance under 
SU(2) conjugation (which rigidly translocates the quadri- 
lateral on the sphere S^) imply in turn the area formula in 
Eq. (|28p , namely that the second argument of T should 
necessarily be proportional to the area of the triangle 
ABC. 

Now we turn to the connection between geometric 
phase for two-level systems and the considerations of 
turns presented above. While geometric phase became 
popular owing to the seminal work of Berry [s^l , it had 
been " anticipated" by Pancharatnam (33| , as pointed out 
by Nityananda and Ramaseshan |34| and subsequently 
by Berry (35| . Some of the other works which could be 
viewed, in retrospect, to have anticipated the geometric 
phase are listed in Rcf. [s^ . 

In the course of his interference experiments with po- 
larized light Pancharatnam faced this question: Given 
two distinct (nonorthogonal) polarization states repre- 
sented by linearly independent Jones vectors E{rii) and 
E{n2), when should one say that these Jones vectors are 
in phase? Motivated by his experiments, Pancharatnam 
arrived at the following answer: E{hi) and E{n2) are in 
phase if and only if the inner product E{ni)^ E{n2) is real 
positive. In other words, being in phase is synonymous 
with maximal constructive interference. He noted that 
"being in phase" defined in this manner is not an equiv- 
alence relation, for it fails the transitivity requirement: 
E{ni) being in phase with E{n2) and E{h2) being in 
phase with E{n-s) docs not necessarily imply that E{n-s) 
is in phase with E{ni). Indeed, he showed that this fail- 
ure in respect of transitivity is geometric in nature, in the 
sense that if E{ni) is in phase with E{h2) and E{h2) is 
in phase with Ein^), then E{h3) will be necessarily out 
of phase with E{ni) precisely by half the area of the 
spherical triangle defined by vertices hi, n2, ^3 G V. 

As a simple illustration of this failure of transitivity, 
consider on the Poincare sphere three points P, Q, R 
identified by unit vectors hp = (0, 1, 0), hq ~ (0, 0, 1), 
and hji = (1, 0, 0) corresponding, respectively, to linear 
polarization at an angle 7r/4 to the xi axis, RCP, and 
linear polarization along the Xi axis. The corresponding 
three Jones vectors may be taken to be 



E{hp) 



1 / 1 



V2 VU' 
E{hp) 



Einq) 



-i ir/4 



1 



1 



(30) 



We have chosen the phase factors multiplying E{hQ), 
E{hji) so as to ensure that E{hp) is in phase with E{hQ) 
and E{hQ) in phase with E{hp). It is readily verified 
that E{hp) is indeed out of phase with E{hp) by 7r/4, 
half the area of the spherical triangle PQR. 

Given any three points hi, h2, h^ on the Poincare 
sphere, we may write out this failure of transitivity in 



the transparent form 

s.Yg[E{hi)^ E{h2) E{h2)^E{h3) E{h3)^E{hi) ] 

= arg [tr ( E{hi)Eihi)^ E{h2)E{h2)^ E{hs)E{h3)^ ) ] 

= arg [tr ( p(ni) ^(na) p(n3 ) ] 



- A3(ni, n2, 1x3). 



(31) 



That A3(ni, h2, h^) defined as above equals the area of 
the spherical triangle with vertices a.t hi, h2, h^ ^ V is 
a result due to Pancharatnam. It may be noted that, 
A3(ni, n2, h^), the argument of the product of pairwise 
inner products of " successive" states, is manifestly gauge 
invariant in the sense that if the Jones vectors E{hj) are 
replaced with e^°'^E{hj), where aj, j = 1,2,3 are arbi- 
trary phases, A3(ni, h2, h^) remains unaffected [every 
vector enters the expression as a bra and as a ket, once 
each ] . It is in view of this gauge invariance that A3 may 
be called a geometric phase. 

If we are given four states corresponding to points 
hi, h2, hs, h4 & V, the associated gauge-invariant ge- 
ometric object is 

arg [E{hi)^Eih2)E{h2)^Eih3)E{h3)^E{hi)Eihi)^E{hi)] 
= axg[E{hi)^ E{h2)E{h2)^ E{h3)E{h3)^ E{hi)] 
+ &rg[E{hi)^E{h3)E{h3)^E{h4)E{h4)^E{hi)] 

= ^A4{hi, h2, hs, hi), (32) 

where in the first step we simply inserted into 
the expression E{h3)'^ E{hi) E{hi)^ E{h3) — 
\ E{h3)^ E{hi)\'^ > 0, which does not affect the 
phase. That A4(ni, 7x2, h^, 714) has to be (a multiple 
of) the area of the quadrilateral follows also from the 
additivity 

Ai{hi,h2,hj„hi) = A3(ni,n2,n3) + A3(ni, n3, n4) (33) 



demonstrated in Eq. ([32]) and independent of Pancharat- 
nam. That our considerations generalize to n states and 
the associated 71-sided polygon is clear. 

In the course of his celebrated proof of the Wigner 
theorem on symmetry in quantum theory Bargmann f36l | 
used the gauge-invariant expression A3 to discriminate 
between unitary and antiunitary symmetries; and it is 
in honour of Bargmann that the authors of Ref. [13, 
named these invariants the Bargmann invariants. In- 
deed, these authors showed that a very general theory of 
geometric phases can be formulated entirely on the basis 
of the Bargmann invariants [13, 113 • The Gouy phase, 
the phase jump a focused light beam suffers at the fo- 
cal point, turns out to be a Bargmann invariant [39| . and 
this could probably be the earliest instance of a geometric 
phase observed in a laboratory. The connection between 
geometric phase and Bargmann invariants has been fur- 
ther explored in Refs. [SmI. 

With this preparation we are now ready to bring out 
the relationship between Hamilton's turns and the ge- 
ometric phase through Bargmann invariants. We be- 
gin with two elementary observations. While "being in 
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, < ip < TT. On the Poincarc sphere this 



phase" is not an equivalence relation in general, on any 
geodesic arc of length < tt on T' it can indeed masquerade 
as one. For proof it suffices to note that the assertion is 
true for the particular case of Jones vectors of the form 
cos((^/2) 
sin((^/2) 

family occupies on the equator all points with azimuthal 
coordinate ip varying from up to (but not including) 
TT. The fact that all these vectors (which correspond to 
linearly polarized states) are in phase with one another 
is obvious. That the claim in respect of masquerading 
applies to an arbitrary geodesic (of extent < tt) on V 
follows from the fact that all such geodesies on V arc 
unitarily equivalent and from the fact that the very no- 
tion of being in phase is unitarily invariant, since it is 
defined through inner products. 

As for the second observation, recall that a turn 
T(n, a) = exp [—ia n ■ r] acts on the Poincarc sphere 
V as rotation of amount 2a about the directed axis n. 
This may be viewed as continuous evolution for a time 
duration a under the constant Hamiltonian h ■ t. Un- 
der this rotation or continuous evolution, states on V are 
driven on circles of constant latitude about n. One of 
these orbits is a great circle. Evolution on these orbits, 
or circles of constant latitude, is not an in-phase evolu- 
tion in general. The geodesic orbit is the only exception: 
A state on V located orthogonal to n evolves in such a 
way that successive states are in phase with one another. 
As an illustration, assume ri = (0, 0, 1) so that 



T(n, Q;)ji=(o,o,i) ^ cxp [-iara] 



cos a — sm a 
sin a cos a 



(34) 



The distinguished great circle in this case is the equator 
and the relevant states are again the linearly polarized 
states considered under the first observation above. That 
the 2x2 real Jones matrix in Eq. p4p drives these real 
Jones vectors in an in-phase manner is obvious. That 
our claims hold for a general n and the associated great 
circle follows from unitary equivalence. This may be 
loosely paraphrased as follows. Under the unitary evo- 
lution driven by r(n, a), states on the great circle or- 
thogonal to n evolve, but not their phases. States on 
the other constant latitude circles evolve with a corre- 
sponding evolution of phases. The distinguished states 
at ±n G 7^ do not evolve; their phases alone evolve by 
±a. 

The connection between turns and Pancharatnam or 
geometric phase emerges quite simply when we combine 
these two observations. Consider the geodesic pentagon 
ABCDE of states shown in Fig.H Let A', B', C, D', E' 
be the mid points of, respectively, AB, BC, CD, DE, and 
EA, and let A be the area of ABODE. Starting with the 
state corresponding to the point A, turn AA' drives it in 
an in-phase manner along AB to B, turn BB' drives B 
in an in-phase manner along BC to C, turn CC drives 
to D along CD, turn DD' drives D to E, and finally 




FIG. 6: Shown, in the case of a geodesic pentagon on V, is the 
connection between Pancharatnam phase and Hamilton's turns. 
A', B', C', D', E' are the mid-points of the sides of the pentagon 
ABODE. On the one hand turns AA', BB', OC, DD', EE' evolves 
or drives the initial state A in an in-phase manner along the five 
sides of the pentagon, resulting in phase change (Pancharatnam 
phase) equaling half the area of ABODE, on the other hand, these 
turns compose to a turn about A of turn length equaling half the 
area of ABODE. 



turn EE' drives E back to A along EA in an in-phase 
manner. 

Now this closed-circuit evolution can be viewed from 
two perspectives. Since the state A is taken along the 
pentagon ABODE back to A in an in-phase manner, 
according to Pancharatnam the final state at A will 
differ from the initial one by a phase equal to A/ 2. 
From the perspective of Hamilton and his turns, we 
have that turn AA', turn BB', turn 00', turn DD', and 
turn EE' acting in that sequence have the combined ef- 
fect of leaving A invariant. That is, A is a fixed point of 
turn AA' -\- turn BB' -\- turn CC -\- turn DD' -\- turn EE'. 
However, any turn leaving A invariant should be a turn 
about fiA, the unit vector specified by As V. Since we 
know from the area formula ([28]) adapted to the present 
case that the turn length of this composite turn is A/ 2, 
we conclude 

T{nA, A/2) = turn AA' + turn BB' + turn CC 

-t-turn DD' + turn EE', (35) 

and this completes the connection between Hamilton's 
turns and Pancharatnam's geometric phase. 

Our demonstration has been for the case of a pentagon, 
but it should be clear that the conclusion generalizes to 
?i-sidcd polygons and, as a suitable limit, to any closed 
circuit on V. 



VI. ELEMENTARY EXERCISES IN THE USE 
OF TURNS 

In this section we illustrate the use of turns in several 
simple situations involving single-qubit gates. The in- 
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sight developed through these elementary exercises will 
prove to be of much value in our analysis to be taken up 
in the next section. 




S 



FIG. 7: The relationship between the Euler parameters C of 
an SU(2) gate n(^, r], Q and the positional spherical coordinates 
on T of the associated turn. The positional coordinates of B, C, A 
arc, respectively, {tt/2, ipi), {0, 'P2), (i" — 2ipi — 1/32). 

Clearly, the first requirement for the effective use of 
turns as a tool-kit for handling SU(2) gates is an abil- 
ity to translate freely between the positional coordinates 
of a turn on the sphere T on the one hand and the 
Euler parameters (angles) of the associated SU(2) ma- 
trix on the other. Given an SU(2) gate u{^,r],(), we 
slide the representative arc of its turn on its great cir- 
cle so that the tail (or head) is on the equator. Let 
the spherical coordinates {6, ip) of the tail and head of 
turn BC be, respectively, {n/2,ipi), {6,ip2), as in Fig.[7l 
It is clear that turn BD = turn EB corresponds to opti- 
cal rotator R{2ip2 ~ 2<^i), turn DC to birefringent plate 
C-T/4+ip2/2('^-26'), and turn AE to C_^/4_|_^j_^2/2(7r - 
26). Thus, we have the suggestive resolution of turn BC 
into its "vertical" (birefringence) and "horizontal" (opti- 
cal rotation) parts : 

ui^,V,0 ~ turn BC = turn BD + turn L»C 

= C_,/4+^,/2(7r-20)i?(2(p2-2v?i). (36) 

In place of this rotation followed by birefringence de- 
composition we could have equally well considered the 
birefringence followed by rotation decomposition. Since 
turn AB = turn BC, the positional coordinates of A are 
(tt — 9, 2ipi ~ ^2)- We have 

V, C) ~ turn AB = turn AE + turn EB 

= R{2^2 - 2^1) C_,/4+^i-^,/2(7r - 20)(37) 

Comparing either of these decompositions with Eq. (|15p 
one readily deduces 

? = -2+^2, 
77 = 11-26, 

C = |-t-<^2-2(^i, (38) 



where ipi,ip2,0 are the positional spherical coordinates of 
the turn associated with m(^, 77, (). 

These are precisely the kind of relationships we were 
after, and it is significant that these expressions which 
connect the positional spherical coordinates on T of a 
turn to its Euler angles are linear. 
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FIG. 8: The commutation relation between a QWP and a HWP, 
yielding the rule for going from the Q-H configuration to the H-Q 
configuration and vice versa. 

Our next exercise concerns the composition of a QWP 
and a HWP, as shown in Fig. [5] Points D, B, E are on 
the equator, A and C are on the circles of 45° latitude, 
and N, S are the polar points. Let <pi, ip2, ^3 be the az- 
imuthal coordinates of the equatorial points D, B, E and 
assume DB = BE or, equivalently, ips = 2ip2 — fi- It is 
clear that turn AS represents Qt/4+i^i/2 while turn NC 
represents (5^/4+1^3/2 • Further, turn SB = turn BN rep- 
resents i?_^/4+,^2/2- 

Now consider the pair of spherical triangles ASB, CNB. 
The angle at N equals the angle at S, in view of the 
assumption ips — f2 = ~ fi- Further, AS = CN and 
SB = NB. Thus, these triangles are congruent, showing 
that turn AB = turn BC. In other words, 

turn AS + turn SB = turn BN + turn NC; 

that is, 

H-Tr/4+V2/2QTr/4+ipi/2 ~ '3ir/4+ip3/2-ff-7r/4+ip2/2 ; 

with ipi + ip3 = 2ip2- (39) 

Denoting tt/A + ipi/2 = ip and — 7r/4 + 1^2/2 = p , the 
constraint tpi + ip^ = 2p2 is equivalent to 7r/4 + (/32/2 = 
2p — p mod TT, and so we have 

H^'Q^ = Q2^'-^H^,, ^p, p'. (40) 

This "commutation relation" shows that the H-Q config- 
uration cannot have a capability not shared by the Q-H 
configuration. 
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FIG. 9: Realization of variable optical rotator using a pair of 
HWPs. This can also be viewed as depicting the special ability 
of a HWP to "absorb" an arbitrary amount of optical rotation. 

The composition of a pair of HWPs is our next ex- 
ercise and this is depicted in Fig. [HI Let (/?2 be 
the azimuthal coordinates of the equatorial points A,B. 
Then turn AB = R{2ip2 - 2(pi), turn SA = turn AN = 
i?_^/4+.^i/2, and turn NB = turn BS = -ff^/4+i/=2/2- One 
readily reads out from Fig.lH] 

turn AB = turn AN + turn NB 
that is, i?(2<^2 - 2(/3i) = i^- 7^/4+^1/2- (41) 

With 7r/4 + (^2/2 = and — 7r/4 + ipi/2 = if we have 

H^H^> ^R{27T + 4{^-^')). (42) 

Recah from Eq. dTS]) that it is i?(47r), and not i?(27r), 
that equals the SU(2) identity tq. [Indeed, R{2n) ~ —tq.] 
Noting that -ff±7r/2+(^ is the inverse of H^, we may rewrite 
the last identity in the form 

i^^^±V2+^' -^(4(^-^5') )■ (43) 

The preceding identity [Eq. (|43)) ] shows that a variable 
optical rotator can be simply realized with a pair of 
HWPs, the effective rotation or optical activity being 
linear in the relative orientation (of the fast axis) of the 
HWPs. 

There is another instructive and important manner in 
which Fig.ini can be read. Since turn BN corresponds 
to i?-7r/4+ip.3/2j the fact that turn AN = turn AB + 
turn BN reads i?-^/4+v3i/2 = ^^-7r/4+<p2/2 R{2(p2 - 2<y9i), 
or HipR{a) ~ Hip-al'i- Similarly, the visual iden- 
tity turn SA + turn AB = turn SB reads R{2(p2 — 
2(^i) i7_^/4+^j/2 = -ff-7r/4+v2/2- We have thus proved 

R{a) H^p = -ff^+Q/4, R{a) — H^_a/i- (44) 

These two identities are equivalent to, and consistent 
with, one another in view of the defining property 



R{a)H^R{-a) = *(a/2) $(-a/2) = i?^+„/2 of 
R{a), and they exhibit the special capability of a HWP 
to 'absorb' optical rotation and still remain a HWP. This 
absorption property, combined with the earlier noted 
property that a pair of HWPs is simply equivalent to 
an optical rotator, implies that three HWPs can fare no 
better than one HWP. 

Our last result shows that any number of HWPs can- 
not, by themselves, realize much portion of the manifold 
of SU(2) gates; indeed, an odd number of them is no bet- 
ter than just one HWP, while an even number is simply 
equivalent to a (variable) optical rotator. In either case, 
not more than a one-parameter family of SU(2) gates gets 
realized. 
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FIG. 10: The subset of unitary gates that can be realized with 
one QWP and one HWP. Shown also is the fact that addition of a 
second HWP can not in any manner enlarge this subset. 

One is thus led to ask how much portion of the SU(2) 
manifold of unitary gates can be realized with two HWPs 
and one QWP. Let us begin with just one HWP and one 
QWP, as shown in Fig.[TUJ where the point A lies on the 
45° latitude circle, and B and C on the equator. Let ipi 
be the azimuthal coordinate of A and C and (p2 that of B. 
Then turn AN corresponds to Q-Tr/4+ipi/2 and turn NB 
to H^/4_^^^/2- As a consequence of the visual identity 
turn AB = turn AN -|- turn NB we deduce that turn AB 
corresponds to if7r/4+v2/2'5-7r/4+v;i/2- It follows that the 
turns or SU(2) gates realizable with one HWP and one 
QWP are distinguished by the property that such turns 
extend from one of the 45° latitude circles to the equator 
or, equivalently, from the equator to a 45° latitude circle. 

However, from the identity turn AB = turn AC -I- 
turn CB we see that such a turn is equivalent to a QWP 
followed by an optical rotator. Since an optical rota- 
tor and a pair of HWPs are equivalent, we conclude 
that turns realizable by two HWPs and one QWP are 
certainly of the same form as turn AB, that is, extend- 
ing between a 45° latitude circle and the equatorial cir- 
cle. Since such turns arc fully parametrized by their az- 
imuthal coordinates on these two circles, we conclude 
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that two HWPs and one QWP can realize only a two- 
parameter family of SU(2) gates, and that this family is 
precisely the one realized using one QWP and just one 
HWP. 

It turns out that with two QWPs and one HWP one 
can realize the entire SU(2) manifold of unitary gates 
[that three QWPs by themselves will not suffice is readily 
seen from the fact that they cannot realize, for instance, 
any SU(2) clement whose turn has length > 3tt/4\. Be- 
fore we turn to the next section for proof of this claim 
of three component realization of all SU(2) gates we ex- 
amine, as our last exercise in this section, the manner 
in which a pair of QWPs transform a variable optical 
rotator into a variable birefringent plate. 
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FIG. 11: Diagram showing that a pair of QWPs can convert vari- 
able optical rotation into variable birefringence. 

The points C, G of Fig.[TT]arc on the 45° latitude cir- 
cles and points A, B arc on the equator. Let AB = 
T]/2. Point K on the geodesic CBG is so chosen that 
CB = KG. It is clear that turn GN corresponds to QWP 
Qo, turn CA to its inverse, Q-j^ji^ and turn AB to R(r\). 

Let us now consider the pair of spherical triangles 
CAB, GNK. We have AC = NG, angle at C = angle 
at G, while CB = KG by construction. Thus, the tri- 
angles are congruent. This means, on the one hand, 
that angle at N = angle at A = 7r/2 so that turn KN 
is a birefringent plate C_7r/4(-) and, on the other hand, 
that KN = AB = 7y/2, so that turn KN corresponds to 
C_^/4(?7). We have thus established 

turn KN = turn KG + turn GN 
= turn CB + turn GN 
= turn CA -I- turn AB + turn GN ; 
that is, C_^/4(r/) = Qo QitI2- (45) 

Conjugating by $(7r/4) we have Co{r]) = 
Q-n/AR{v)Q~'7T/ii which on conjugation by ^{tp) 



yields 

(46) 

This relationship between variable optical rotation and 
variable birefringence is what we set out to demonstrate. 
However, to the extent that the variable optical rotator 
is not a preferred component, the fact remains that this 
may not be the most convenient way to realize variable 
birefringence. 

VII. REALIZATION OF ALL SU(2) GATES 
USING TWO QWPS AND ONE HWP 

Since HWP and QWP have each only one (rotational) 
degree of freedom, and since SU(2) is a 3-parameter man- 
ifold, it is clear that at least three components are re- 
quired to realize even a small nontrivial (nonzero mea- 
sure) part of this manifold. We have already noted that 
one QWP and two HWPs cannot fare any better than a 
QWP plus a HWP. In the present section we prove, us- 
ing insights developed through the elementary exercises 
of the last section, that two QWPs and one HWP are suf- 
ficient to realize the full group manifold of SU(2) gates. 

The proof is straight-forward and relies entirely on 
Fig. [121 The points A, B, D are on the equator, N and S 
are the polar points, and C, G are on the circles of 45° 
latitude. The point F, the intersection of line (great circle 
arc) DN with line CBG, is not assumed to be on the circle 
of 45° latitude; the fact that F indeed lies on this circle 
will emerge by itself. The point K on the line CBG is so 
chosen that turn CB = turn KG. We do not assume that 
the angle GNK equals 7r/2. The fact that NK and NG 
are orthogonal at N will unfold on its own; turn KN will 
then correspond to a birefringent plate C_^/4(?7 ), where 

T] equals twice the arclength of KN. It is this turn cor- 
responding to variable birefringence that will eventually 
become the focus of our attention. 

We assume AB = BD = rj/2. Our first task is to prove 
that the point F lies on the 45° latitude circle. We begin 
by noting that turn GN corresponds to QWP Qo and so 
is also turn CS. 

Consider the pair of spherical triangles ABC, DBF. 
The angle at A equals the angle at D (both equal 7r/2). 
Both triangles have the same angle at B, and AB = BD 
by construction. Thus, the two triangles are congruent. 
It follows that DF= AC = 7r/4, showing that F lies indeed 
on the 45° latitude circle, and that turn BF — turn CB. 
Since F is on the 45° latitude circle, turn NF corresponds 
to QWP Q,,/2. Since AB = 77/2, turn BN = turn SB 
corresponds to HWP H^/2+,^/i = -ff-7T/2+r;/4- 

Now consider the pair of spherical triangles ABC, 
KNG. The angle at C obviously equals the angle at 
G. Further, CB = KG by construction and AC = GN 
(= 7r/4), showing that the two triangles are congruent. 
As one consequence we see that the angle at N is 7r/2, 
showing that turn NK is indeed the birefringent plate 
C'-ir/4( ■ )• another consequence we have KN = AB = 
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FIG. 12: Realization of all SU(2) gates using just two QWPs and 
one HWP. Two perspectives are presented for visual convenience. 

BD = 77/2, showing that turn KN = C_7r/4('?)- But, in a 
visuahy obvious manner, 

turn KN = turn KG + turn GN 
= turn BF + turn GN 
= turn BN + turn NF + turn GN. (47) 

Since turn KN = C_^/4(77), and since the three turns 
on the right-hand side of the last equation equal, respec- 
tively, iy7r/2+r)/4i '3?7/2: ^-nd Qq, we have proved 

C-TT/iiv) = Qo Qr,/2 HT,/2+rj/4, (48) 

demonstrating the realizability of variable birefringence. 
We see from Fig. [12] that turn KN could have also been 
developed in a somewhat different but equivalent man- 
ner: 

turn KN = turn KG + turn GN 
= turn CB + turn GN 
= turn CS + turn SB + turn GN. (49) 



It is clear that this corresponds to the multiplicative iden- 
tity 

C-tt/4{v) = Qq H^i2+n/A Qo- (50) 

Incidentally, the fact that these two realizations of 
C-Tr/4{v) respectively in the Q-Q-H and Q-H-Q configu- 
rations are equivalent can also be verified using the H-Q 
commutation relation in Eq. (|40p . 

Conjugating Eq. (|48| by $(7r/4), which corresponds to 
rigidly rotating Fig. [T^ by 7r/2 about the polar axis, we 
have 

Coiv) = Qtt/4 <37r/4+r,/2 -ff-7r/4+J7/4, (51) 

a variable bircfringcnt plate with the fast polarization 
along the xi axis. The reader will appreciate that Fig.fT^ 
was crafted for C_tj/4^{vi) rather than Co{r]) simply for 
visual clarity. 

Conjugating the last equation by $(^/2), which 
amounts to a further rigid rotation of Fig.[T2] by ^ about 
the polar axis, one obtains 

/4+e/2 Q7r/4+4/2+r;/2 -ff-7r/4+J/2-|-r,/4 • 

(52) 

Right multiplying by the optical rotator + (), and 
using the special ability of HWP to "absorb" rota- 
tion, namely H^.)R{a) ~ i?(.)_Q/4, we have in view of 
<S., V, C) = Q/2(??)i?(e + C) given in Eq. m 

^{^7 Vi C) = 07r/4+e/2 Qir/4+C/2+r,/2 ^-7r/4+(e+i)-C)/4, 

(53) 

showing that all SU(2) gates can be realized using 
two QWPs and one HWP in the Q-Q-H configuration. 
Uniqueness of this realization as well as the fact that the 
required (angular) positions of the plates are linear in the 
Euler angles should be appreciated. 

That a similar assertion holds for the other two pos- 
sible configurations, namely, Q-H-Q and H-Q-Q, fol- 
lows immediately from the Q-H commutation relation 
Qo.Hp^HpQ2p-^ of Eq.JlOl): 

C) = Q7r/4-|-5/2 -f^-7i-/4+(4-|-r;-C)/4 Q7r/4-C/2; (54) 
?7, C) = ^-7r/4-(-(5+i)-C)/4 '37r/4+(i)-C)/2 '37r/4-C/2- 

(55) 

We have thus completed a proof of the main result of this 
section, indeed this paper: 

Theorem: All SU(2) gates m(^, rj, C,) can be realized 
with just two QWPs and one HWP equally well in any 
one of the three conceivable configurations Q-Q-H, Q-H- 
Q, or H-Q-Q, as detailed, respectively, in Eqs. (|53)) . (|54l) . 
and (j55p . The realization is unique in each case, and in 
each configuration the angular positions of the plates are 
linear in the Euler angles ^, r;, C, of the gate. 
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FIG. 13: The assembly of the universal gadget which realizes all 
single-qubit unitary gates. The three wave plates Q, H, Q' are 
coaxialy mounted, with a provision being made to assign and read 
the orientations of their fast axes any triplet of values. 



VIII. UNIVERSAL SINGLE-QUBIT UNITARY 
GATE 

The theorem proved above enables assembUng of a 
universal single-qubit unitary gate as described below. 
While the universal gate can be assembled in any one 
of the three configurations Q-Q-H, Q-H-Q, or H-Q-Q we 
choose the configuration Q-H-Q simply in order to be 
concrete. 

Let two QWPs Q, Q' and a HWP H be coaxially 
mounted in a cylindrical shown in Fig.[T31 As- 

sume that a provision is made to rotate each one of the 
three plates Q, H, Q independently about the axis of 
the cylinder, and that provision is made to read out the 
angular coordinates (/3i, Lp2, (^33 of their fast axes on their 
respective (semicircular) dials. We assume further that 
the assembly is so used that light passes through Q, H, 
Q in that order. 

Clearly, the SU(2) matrix corresponding to the entire 
assembly is Q^^ H^^ Qtpi- Thus, to realize any unitary 
gate 77, C) € SU(2) we have to simply arrange these 
dial positions to read 

ipi = 7r/4 — C/2 mod TT, 

V2 = -7r/4+(^ + r/-C)/4 mod vr, 

^3 = 7r/4 + C/2 modTT, (56) 

as may be seen from Eq. ([51]) . 

That the entire SU(2) manifold of unitary gates can 
be realized using a single assembly has its advantages. 
For instance, if the only imperfections of the wave plates 
from ideal ones are residual losses, the fact that the 
total loss of the assembly is the same for realization 
of every SU(2) gate, independent of the Euler angles 
(^, 77, C,) of the gate, implies that it can be accounted 
for more easily. Further, the fact that the dynamical 
phase through the system remains the same for all gates 
can prove to be of particular importance in interference 



considerations, particularly in the context of geometric 
phase experiments [3, • 



IX. CONCLUDING REMARKS 

In this paper we have developed in considerable detail 
the pictorial construction of Hamilton into a tool-kit for 
handling problems of synthesis and analysis in situations 
where the unitary group SU(2) plays a central role. This 
group pervades nearly all areas of science, either directly 
or through its cousin S0(3). Thus, this tool-kit should 
be of interest to a wider audience beyond quantum com- 
putation and quantum information. [In particular, the 
formalism and results presented here should be of much 
interest to classical polarization optics.] It is for this rea- 
son that we have developed Hamilton's construction itself 
in sufficient detail, taking care to bring out its connection 
with geometric phase and the Bargmann invariants. It is 
in view of this detailed preparation that we believe the 
manipulations with turns presented in Sections 6, 7, and 
8 will be found to be fully accessible to a broad spectrum 
of readership. 

As noted in the Introduction, the central result pre- 
sented as a theorem towards the end of the paper is not 
new in itself. Our presentation is fashioned to work to- 
ward this result for two reasons. First, to demonstrate 
how this geometric approach is suggestive and visually 
transparent compared to the algebraic approach. Sec- 
ond, this result acts as a focal point in the sense that 
in working towards it most of the simple manipulations 
with turns are conveniently developed and demonstrated 
in stages. 

It is true that all the results developed here geometri- 
cally could be algebraically verified through matrix mul- 
tiplication. However, it is in the geometric representation 
that synthesis results suggest themselves in a vivid or vi- 
sual manner. We may cite the H-Q commutation relation 
in Fig. 8 and the "absorption" of optical rotation by a 
HWP (as also the realization of variable optical rotator 
using a pair of HWPs) in Fig. 9 as effective illustrations 
of this advantage. 

Finally, our presentation of turns in this paper is in 
the context of the unitary group SU(2). The reader will 
easily realize that this geometric representation readily 
translates to the rotation group S0(3) if a turn of length 
I and the reversed turn of length tt — £ are identified. 
This amounts to identifying the null turn with the turn 
of length TT, which is clearly the same thing as identifying 
U of SU(2) with - [/ of SU(2). Thus, the turn length ^ in 
the case of S0(3) gets restricted to the range < ^ < 7r/2. 
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